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TECHNICAL  REPORT 


The  main,  long  term,  objective  of  this  work  is  research  and  development  of  the 
theory  and  constructive  computational  algorithms  for  synthesis  of  single  and  dual 
reflector  anteima  systems.  The  work  is  based  on  analytic  and  numerical  procedures  for 
solving  the  underlying  nonlinear  boundary  value  problems. 

To  this  end  we  have  been  carrying  out  investigations  in  the  following  directions: 

1.  Conditions  for  solvability  of  the  direct  problem  of  reconstructing  the  reflector  antenna 
with  uniform  output  density  of  the  reflected  rays  from  the  given  far  field,  input  aperture, 
and  prescribed  in  advance  (non-radially  symmetric)  input  density. 

2.  Properties  of  reciprocal  reflector  and  connections  between  the  direct  and  inverse 
problems. 

3.  Construction  and  testing  of  an  algorithm  based  on  a  "diffusion"  -  type  scheme  for 
solving  numerically  the  boundary  value  problem  formulated  in  1. 

4.  Formulation  of  the  dual  reflector  problem  as  direct  and  inverse  boundary  value 
problems  and  investigation  of  appropriate  solvability  conditions. 

In  item  1  we  succeeded  in  deriving  the  equation  in  concise  and  explicit  form.  As  far  as  we 
know,  in  such  form  this  equation  is  obtained  for  the  first  time.  All  the  terms  in  the 
equation  have  a  simple  and  clear  geometric  meaning  and  can  be  computed  numerically 
by  efficient  numerical  procedures.  Our  derivation  is  based  on  a  general  procedure 
utilizing  differential  geometric  methods.  We  also  use  the  same  methods  for  deriving  the 
equations  of  the  dual  reflector  problem  in  item  4.  We  have  also  shown,  in  explicit  form, 
the  connection  between  the  inverse  and  direct  problems  in  the  single  reflector  setting 
(item  2).  Using  this  connection  we  established  conditions  for  solvability  of  direct 
problem.  Regarding  item  3  we  developed  and  tested  an  algorithm  for  solving  numerically 
the  direct  single  reflector  problem.  Our  approach  here  is  via  a  certain  special  "diffusion" 
-type  procedure.  With  this  approach  it  is  possible  to  avoid  a  costly  numerical  inversion  of 
a  highly  nonlinear  second  order  differential  operator. 

The  results  are  being  organized  in  a  series  of  papers.  Two  of  them  "Differential- 
geometric  methods  in  design  of  single  and  dual  reflector  antennas"  and  "On  one  direct 
problem  in  the  reflector  antenna  theory"  are  complete  and  submitted  for  publication. 
Copies  of  both  papers  are  attached  to  this  report.  The  third  paper  "On  the  theory  of 
dual  offset  reflector  antennas"  is  being  prepared  for  publication. 


Differential-geometric  methods  in  design  of  single  and  dual 

reflector  antennas* 

by 

Vladimir  Oliker,  Elsa  Newman®  ,  and  Laird  Prussner 

Introduction 

This  paper  is  the  first  in  a  series  of  three  papers  in  which  we  study 
the  theory  and  numerical  methods  in  synthesis  of  reflector  antennas. 

The  problem  of  synthesizing  a  single  or  dual  reflector  antenna 
sgstem  producing  a  pre-specified  intensity  distribution  on  a  far-field  or 
on  a  target  domain  continues  to  attract  considerable  attention  of 
researchers  and  practitioners.  In  the  geometric  optics  approximation 
(GO)  the  basic  laws  of  reflection  can  be  used  to  derive  a  system  of 
three  first  order  partial  differential  equations  (PDE's)  corresponding  to 
the  problem.  This  was  done  in  the  60s  by  B.  Kinber  [7]  and  V.  Galindo- 
Israel  [4].  This  system  of  PDE's,  roughly  speaking,  consists  of  two 
"parts";  two  equations  in  the  system  express  the  Snell  law,  and  the  third 
equation  relates  the  output  intensity  distribution  to  the  intensity  of  the 
primary  source.  The  latter  is  a  strongly  nonlinear  equation  which  is 
essentially  a  condition  on  the  Jacobian  of  the  "reflector"  map 
transforming  the  input  spherical  wave  front  into  the  output  front. 

It  is  well  known  (see,  for  example,  {3],  Ch.  I)  that  a  first  order 
system  of  PDE's  admits  a  smooth  solution  only  if  certain  integrability 
conditions  are  satisfied.  It  was  observed  by  Kinber  [7]  and  Galindo- 
Israel  and  Mittra  (see  [5]  and  further  references  there)  that  for  the  first 
order  PDE's  system  describing  reflector  antennas  these  integrability 
conditions,  in  general,  may  not  be  fulfilled.  In  fact,  it  is  not  difficult  to 
write  down  explicitly  the  integrability  conditions  for  the  two  equations 


'This  research  v/as  supported  by  AFOSR  under  contract  F499620-91- 
C-0001.  The  United  States  Government  is  authorized  to  reproduce  and 
distribute  reprints  for  governmental  purposes  notwithstanding  any 
copyright  notation  herein. 
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expressing  Snell's  law;  see  [5],  [61.  However,  this  is  not  sufficient,  since 
the  third  equation,  relating  the  input  and  output  intensities,  is  not  taken 
into  account.  Consequently,  one  can  not  expect,  in  general,  that  the 
problem  will  have  a  solution. 

The  exceptional  situation  when  the  integrability  conditions  for  the 
entire  system  are  always  satisfied  is  the  case  of  axially  symmetric 
reflectors  and  axially  symmetric  data.  It  has  been  known  for  some  time 
and  was  shown  rigorously  in  [15]  that  in  this  case  one  can  give  explicit 
necessary  and  sufficient  conditions  for  solvability  of  the  problem.  In  the 
subsequent  paper  we  will  show  that  if  the  input  domain  and  the  far- 
field  are  coaxial  and  there  is  an  axially  symmetric  solution  of  the 
problem  with  axially  symmetric  intensities  then  one  can  always  find  a 
solution  of  the  first  order  PDE's  system  for  prescribed  non-axially 
symmetric  intensities  close  in  some  norm  to  the  axially  symmetric  one. 

In  reflector  antenna  systems  with  two  or  more  reflecting  surfaces 
there  are  additional  degrees  of  freedom  to  control  the  output  amplitude 
and  phase.  However,  the  more  reflecting  surfaces  are  in  a  system  the 
larger  the  number  of  PDE's  describing  it.  Consequently,  the  complexity  of 
the  system  increases  and  its  analysis  becomes  harder.  As  far  as  we 
know  up  to  present  there  are  no  rigorous  results  concerning  such 
systems.  The  authors  usually  use  numerical  computations  to  obtain  some 
acceptable  variant  of  a  "solution"  (even  when  a  true  solution  may  not 
exist).  In  the  most  recent  v/orks  of  K.  Shogen,  R.  Mittra,  V.  Galindo- 
Israel  and  W.A.  Imbriale  various  difficulties  in  numerical  treatment  of 
these  equations  are  reported  [17]. 

In  1972  J.  Schruben  [16]  considered  the  problem  of  designing  a 
lighting  fixture  which  would  produce  a  pre-specified  intensity  pattern 
on  a  plane  aperture.  In  her  approach,  the  reflector  surface  is  described 
as  a  graph  of  a  scalar  function  over  a  domain  on  a  unit  sphere  (the  input 
aperture)  and  the  Jacobian  of  the  reflector  map  is  expressed  in  terms  of 
this  function  and  its  derivatives  up  to  second  order.  The  corresponding 
expression  is  a  second  order  nonlinear  PDE  of  Monge-Ampere  type. 
Unfortunately,  the  use  of  a  specific  coordinate  system  made  the 
formulas  in  [16]  quite  complicated,  and  Schruben  does  not  even  derive 
explicitly  the  equation  of  the  problem.  Her  main  concern  was  to 
describe  conditions  for  ellipticity  of  the  operator  involved. 

The  principal  advantage  in  this  approach  is  that  instead  of  a  system  one 
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has  to  deal  here  with  one  equation  for  one  scalar  function.  This 
equation,  however,  is  of  second  order  and  strongly  nonlinear. 

In  1974  B.  Westcott  et  al.  initiated  a  different  approach  to 
synthesis  of  reflectors.  In  case  of  a  single  reflector,  required  to  produce 
a  prescribed  far-field  pattern,  this  approach  to  the  problem  corresponds 
to  formulating  it  as  an  inverse  problem.  It  is  based  on  the  observation 
that  the  reflector  can  be  parametrized  in  a  special  way  by  points  on  the 
far-field  (or  target)  and  actually  recovered  from  one  scalar  function. 

The  latter  must  also  satisfy  a  second  order  PDE  of  Monge-Ampere  type 
relating  the  intensity  of  the  primary  source  and  of  the  desired  output. 
This  approach  was  pursued  by  Westcott  and  his  colleagues  in  a  number 
of  publications;  see  111],  11],  [2],  Il9],  [18],  and  further  references  there. 

In  all  of  these  approaches,  rigorous  general  results  concerning 
existence  and  uniqueness  of  solutions  are  lacking.  Eventually,  the 
authors  always  have  to  resort  to  numerical  calculations.  Some  progress 
towards  establishing  rigorous  uniqueness  results  was  made  by  Harder 
[10],  who  considered  the  single  reflector  problem  in  the  setting  of  the 
approach  by  Westcott  et  al.  Also,  in  the  inverse  problem  setting 
rigorous  results  establishing  existence  of  non-axially  symmetric 
solutions  with  prescribed  non-axially  symmetric  densities,  close  (in  a 
certain  norm)  to  axially  symmetric,  were  obatained  by  Oliker  [13]  (see 
also  [14]  for  related  results)  in  the  case  when  the  problem  is  treated. 

In  this  paper,  we  study  the  “direct"  problem  of  synthesizing  the 
reflector  surface.  Our  starting  point  is  similar  to  that  of  Schruben  in 
the  sense  that  we  also  describe  the  sought  reflector  surface  as  a  graph 
over  a  spherical  domain  (the  input  aperture)  and  look  for  a  second  order 
PDE  which  the  function  describing  the  graph  satisfies.  However,  in 
contrast  with  the  paper  of  Schruben,  we  succeed  in  deriving  such  an 
equation  in  concise  and  explicit  form.  This  is  done  in  both  cases  of 
single  and  dual  reflector  antennas.  In  this  paper  we  present  our  results 
for  a  single  reflector  antenna.  The  results  for  dual  reflector  antennas 
will  be  presented  in  a  subsequent  publication. 

The  corresponding  expressions  are  relatively  simple  and,  most 
importantly,  contain  familiar  geometric  quantities  for  which  various 
stable  procedures  for  numerical  computation  are  available.  Our 
derivations  are  based  on  differential  geometric  methods.  We  also  show 
that  in  the  case  of  a  single  reflector  surface  our  equation  and  the  one 
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derived  by  Brickell-Marder-Westcott  [2]  are  connected  by  a  simple 
transformation.  In  fact,  we  show  that  the  reflector  surfaces  described 
in  the  direct  problem  and  the  reflector  surfaces  constructed  in  the 
inverse  problem  are  actually  what  is  known  as  reciprocal  reflectors.  In 
the  complex  analytic  formulation,  such  a  connection  between  the  direct 
and  inverse  problems  was  established  in  [2].  Our  approach  to 
establishing  this  connection  does  not  involve  complex  analysis.  In  the 
real-valued  form  this  connection  becomes  quite  transparent  and  follows 
easily  from  the  formulas  we  develop.  We  also  discuss  the  question  of 
optimal  boundary  conditions  to  be  imposed  on  the  solution.  This 
question  is  important.  As  we  pointed  out  earlier,  the  problem,  in 
general,  may  be  overdetermined  and  therefore  lack  a  solution. 

The  paper  is  organized  as  follows.  In  section  1  we  review  some 
facts  from  elementary  differential  geometry  and  develop  basic 
geometric  formulas  for  reflecting  surfaces.  In  section  2  we  compute  the 
Jacobian  of  the  "reflector"  map  and  relate  it  to  the  input  and  output 
power  densities.  The  main  result  here  is  the  formula  (2.12)  expressing 
the  Jacobian  in  geometric  quantities.  As  far  as  we  know,  in  such 
explicit  form  this  expression  appears  for  the  first  time.  This  expression 
turned  out  to  be  very  useful.  In  section  3  we  use  it  to  construct  the 
reciprocal  reflector  and  show  the  connection  between  the  direct  and 
inverse  formulations  of  the  problem.  In  section  4  we  present  and 
analyze  in  our  framework  the  formulations  of  the  problem  as  a  first 
order  system  of  PDE's  and  as  a  boundary  value  problem  for  a  second 
order  nonlinear  PDE  of  Monge-Ampere  type.  In  section  5  we  give  a 
rigorous  treatment  of  the  axially  symmetric  case  of  the  direct  problem. 
In  section  6  a  numerical  algorithm  based  on  formula  (2.12)  is  presented. 

1.  Preliminaries 

1.1.  In  three  dimensional  space  let  S  be  a  unit  sphere  centered 
at  some  point  0.  Fix  a  Cartesian  coordinate  system  with  the  origin  0. 

Let  be  a  domain  on  S,  Q  its  closure,  and  m  a  unit  vector  with  endpoint 
in  ?2.  Let  p  be  a  positive  function  in  9.  and  set  r(m)  =  p(m)m.  Then  r 
defines  a  surface  F  projecting  radially  from  0  univalently  onto  9. 

Denote  by  n  the  unit  normal  vector  field  on  F  and  assume  that  F  is 
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oriented  so  that  <nn,n>  >  0  everywhere  on  F.  An  illustration  of  this 
situation  is  given  on  the  figure  on  p.  8. 

We  now  recall  some  basic  notions  from  Differential  Geometry  that 
will  be  needed  later.  These  can  be  found,  for  example,  in  E.  Kreyszig  19]. 

Let  (u^u^)  be  some  smooth  local  coordinates  on  S.  Put  u  =  (u^u^). 
Then  m  =  m(u)  (sm(u^u^))  is  a  vector  valued  function  of  u  giving  the 
position  vector  of  any  point  in  c  S.  For  that  reason  m(u)  is  viewed  as 

a  unit  vector  in  and  also  as  a  point  in  Q.  As  usual,  we  put  f(m(u))  s 
f(u’,u^)  for  any  function  f:S“>R.  Everywhere  in  the  paper,  we  use  the 
following  range  of  indices:  l<i,j,k,...<2. 

The  first  fundamental  form  e  =  ejjdu’duJ  of  S  has  coefficients 


where  <  ,  >  denotes  the  scalar  product  in  R^  and  =  9m/9u^  Here  and 

everywhere  below  the  summation  convention  over  repeated  lower  and 
upper  indices  is  in  effect.  The  matrix  [e^j]  is  symmetric  and  invertible; 

its  inverse  is  denoted  by  [e’Jj. 

We  assume  that  the  coordinates  u^  u^  are  chosen  so  that 
<m,nn  jX  m2>  >  0  in  ?2. 

The  area  element  da  of  S  (in  n)  is  given  by 

da  =  I  mi  X  m2|du’du2  =  (m,mi  x  m2)du’du2  =  ydet(e  ..) 

The  first  fundamental  form  g  =  g  ydu’duJ  of  the  reflector  surface  F 
has  coefficients 


Qij  =  <n,rj> 


PiPj 


where  =  9r/9u^  and  =  9p/9u\  Additionally,  we  set 


Vp  =  p^g’^rj,  Vp  =  piO^jmj 
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Where  [g’j]  =  [gjjl  \  Then 

|Vp|2  =  <vp,Vp>  =  pjpjg’j,  |Vp!2  =  <Vp,Vp>  =  p^pje’j. 

The  vectors  m^(u),  m2(u),  and  m(u)  form  a  basis  of  at  every 

u  e  Q  and  any  vector  In  can  be  expressed  in  terms  of  this  "moving” 
basis.  In  particular,  we  may  express  the  unit  normal  vector  field  n(u) 
on  the  reflector  surface  F  as 


II  -  — y  -  V  I  .  I  / 

Vp^ *  I'^pI^ 

Since  Tj  =  pjm  +  pm^  we  have  <rj,n>  =  0.  Obviously,  |n|  =  1  and 

so  n  is  indeed  the  unit  normal  field  on  F. 

It  follows  from  (1.1)  that 


<r,n)  =  p{m,n>  = 


(1.2) 


Since  by  our  assumption  <m,n>  >  0  and  ]  r  |  =  p  >  0  on  F,  we  see  that 


<r,n>  >  0  on  F. 


(1.3) 


The  covariant  derivatives  relative  to  the  first  fundamental  form  of 
F  are  defined  for  any  scalar  function  f:  F->R  as  follows: 


V  f  =  f  -  r'^f  , 
ij  ij  ij  k 


where 
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f  =  ^  f  = 

'  9u'  ’  'j  3u'3u^  ’ 

and  the  T's  are  the  Christoffel  symbols  associated  with  g. 

Similarly,  for  Q  c  S  and  f:  Q  ->  R,  the  covariant  differentiation  is 
defined  as 


Vif  -fi, 


where  the  f's  represents  the  Christoffel  symbols  associated  with  the 
first  fundamental  form  e  of  S. 

If  V  =  (v,,V2,V3)  is  a  vector  field  on  F  then 


V^V  =  (VjV,,V^V2,V^V3) 


and  similarly  one  defines  Vjjv,  and  V^v,  if  v  is  a  vector  field  on 

S7  c  S.  The  coefficients  of  thr  second  fundamental  form  b  =  b  jjdu’duJ  on 
the  surface  F  are  given  by 

bij  =  <r^j,n>  =  -  <rj,nj> 

where  r^j  =  32r/3u’9uj. 

According  to  the  classical  derivation  formulas 

r.  =  +  b..n  (1.4) 

ij  ij  k  ij 

1.2  Suppose  a  light  ray  is  originated  at  0  in  the  direction  m  and  is 
reflected  at  the  point  r(m)  in  the  direction  y. 
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By  Sneirs  law, 


y  =  m  -  2  <m,n>n, 

and  we  may  consider  the  map  Q  —>  S,  ?J(m)  =  y(m),  m  e  ?2.  Put 

^  =  tJ(sT).  The  may  transforms  Q  — >  Consequently,  we  can  relate 
the  area  elements  in  H  and  in  its  image  w  by  computing  the  Jacobian  of 
the  map  7$.  Note  that  since  y(m)  is  the  unit  vector  in  the  reflected 
direction,  the  Jacobian  J(m(u))  is  given  by 


(u)\ 

— L  (1.5) 

.(u)> 


J(m(u))  =  ± 


d^^(m(u)) 


da(u) 


(u),y^ 


m.(u),m 


We  assign  a  ±  sign  to  the  Jacobian  according  to  whether  7$  preserves  the 
orientation  or  reverses  it.  This  is  equivalent  to  considering  the 
volumes'  relationship 


J(m)<m,(m  ,xm2)>  =  <y,(yixy2)>  . 
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Since  <m,  nrii  x  m2>  >  0,  J  is  positive  if  <y,  x  y2>  >  0  and  J  is 
negative  otherwise.  Note  that  |  <y,(yixg2)>  |  =  |  det<yj,yj>  | . 


1.3  We  will  need  several  different  expressions  for  the  vector 
functions  m(u)  and  y(u).  First,  ve  find  an  expression  for  m(u)  in  terms 
of  the  basis  r^Cu),  r2(u),  n(u).  These  vectors  indeed  form  a  basis  since 

|ri(u)  X  r2(u)|  =  p2(u){|  vp(u)|2  +  p2(u)}det(e^j(u))  >  0, 

and  n(u)  is  perpendicular  to  ri(u)  and  r2(u). 

Now,  expressing  m(u)  in  terms  of  r,(u),  r2(u),  and  n(u),  we  find 
(omitting  the  argument  u) 


m 


P|g'V  *  7'  -  Ivppn 


(1.6) 


Obviously,  this  formula  is  valid  only  if  jvpl^  <  1.  However,  it  follows 
from  (1.2)  and  (1.6)  that 


<m,n> 


(1.7) 


Since  it  is  always  assumed  that  <m,n>  >  0  and  p  >  0  on  F,  the  condition 
I  Vp  1 2  <  1  is  fulfilled. 

We  may  express  y  in  terms  of  p,  m,  and  their  derivatives.  It 
follows  from  Snell's  law  and  (1.1)  that 


y  =  m 


-  2 


p(pm  -  Vp) 
|7p|2  ♦  p2 


(1.8) 


Using  Snell's  law  and  (1.7),  v/e  may  express  y  without  explicit  use  of  m: 

y  =  - 7'  -  KpI^ 


n 


(1.9) 
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2.  Computation  of  the  Jacobian  of  the  map  7! 
and  the  “balance**  equation 

2.1. Here  we  find  an  explicit  expression  for  the  Jacobian  J(m).  As 
It  follows  from  (1.5),  we  need  an  expression  for  det<y^,yj>. 


Proposition  2.1.  Let  F  be  a  reflector  surface  as  In  section  1.  Put 
H..  =  V..p  +  J\  -  I  Vp|2  b..  (2.1) 

IJ  ^  I  r  I  ,j 

P  -  9  Pk  • 


Then 


det  H. 


IJ 


/det<y.,y  >  =  — r  —  • 

V  '  i  M  9,:  -  |Vp|2 


V 


(2.2) 


Remark.  Note  that  since  |  Vp  |  ^  <  i,  (2.2)  Is  well  defined. 
Proof.  We  begin  by  showing  that 


=yi.y,> 


=  H.  H. 

Ik 


.ks 


k  5 

p  p 


I  -  |Vf.| 


(2.3) 


Using  (1.9),  we  differentiate  y  covarlantly  relative  to  the  form  g  and 
obtain 


..  g  V,.pp.  _ 

Vy=y  =  V  pgV  +  -  .  '= J=  n  -  J]  - 1  Vp  |  ^  n  . 

i  *  j  yi  -  |vp|2  V  I  n  i 


To  simplify  this  expression  we  use  (1.4)  and  the  Weingarten  equations 
[9],  p.  126, 
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"i  =  -  . 

Now  we  see  that 


Calculating  <yj,yj>,  we  immediately  arrive  at  (2.3). 
On  the  other  hand,  we  have 


p'p^  \ 

1  -  |Vp|2/ 


_ 1 _ 

det(g..)(l  -  |Vp|2) 


(2.5) 


Consequently,  evaluating  the  determinants  on  both  sides  of  (2.3),  we 
obtain  (2.2).  The  proposition  is  proved. 

2.2.  It  will  be  useful  to  find  an  expression  of  (2.2)  in  terms  of  p 
and  its  derivatives  on  S. 

Proposition  2.2.  Let  H^j,  b^j,  and  e^j  be  as  before.  Then 
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=  2b^j<m,n>  +  pejj 


(2.6) 


Proof.  It  follows  directly  from  (1.6)  that 

H^j  =  2b|j<m,n>  +  V^jp  (2.7) 

Since  p^  =  <r,r>,  and  pp^  =  <r,r^>,  we  differentiate  covariantly  and 
obtain 

p,Pj  »  pVjjp  =  <r,,rj>  <r.Vi^r>  =  g,j  *  <r,n>. 

Hence,  pVjjp  =  b,j<r,n>  -  p,p(  <•  gji.  But 

9|J  =  *  pfpPjm  ♦  P^j^  =  PiPj  *  P^etj- 


Thus, 

pV|jp  =  b^j<r,n>  +  p^e^j. 


Dividing  by  p  and  substituting  in  (2.7),  we  obtain  (2.6).  The  proposition 
is  proved. 

The  following  expression  for  the  second  fundamental  form  in  terms 
of  p  is  derived  in  [12] 


b.. 

ij 


pV,]P  -  p^e.j  -  2p,p. 


(2.8) 


Combining  it  with  (2.6)  and  (1.2),  we  get  the  following  expression 


(2.9) 
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2.3.  We  now  express  det(gjj)  in  terms  of  det(ejj).  Recalling  that 
g^j  =  <r^,rj>  =  pjpj  +  p^ejj  and  evaluating  the  determinant  of  g,j,  we 
find 


detg^j  =  p2(p2  +  I  vp  1 2)det(e|j). 


(2.10) 


We  use  equations  (2.6),  (2.10),  and  (1.7)  to  rewrite  the  equation  (2.2)  as 


I  det(2b  <m,n>  +  pe  ) 
J_  J _ ij  ij 

p^  ^det(e .) 


(2.11) 


2.4.  Summarizing  the  preceding  results,  we  may  now  obtain  the 
expression  for  the  Jacobian  of  the  map  '6:  Q  ->  Zo.  Namely,  taking  into 
account  the  sign  convention,  we  obtain  from  (1.5)  and  (2.11) 


ydet<y.,y.>  ^  det(2b.  .<m,n>  +  pe..) 


/det<m.,m.>  p^  det(e. .) 

V  '  J  '1 


=  G(p). 


(2.12) 


We  use  expressions  (2.8)  and  (2.H)  to  find 


J  = 


(p2  ♦  |9pr) 


det{2p9,,p  -  (p2  -  |^^pr)e..  -  4p.p.) 

det(e..) 

u 


(2.13) 


2.5.  Using  the  expression  for  J(m),  we  are  now  in  a  position  to 
relate  the  energy  of  the  input  primary  source  emitting  a  power  density 
Km)  to  the  desired  output  povver  pattern  V(y).  Namely,  if  |  dtf(m)  |  is 
the  area  element  in  co  expressed  via  the  map  then  we  have  the  point- 
wise  balance  equation 


V(y(m))  |d?J(m)|  =  I(m)dcJ(m),  m  e 


(2.14) 
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or  by  (1.5) 


V(y(m))  I  J(m)  I  =  I(m),  m  6  ft. 


(2.15) 


Taking  into  account  the  sign  of  the  Jacobian  we  obtain  from  (2.15)  and 

(2.12) 

(b.  .(m)<m,n(m)> 

2-^ - - +  e..(m) 

p(m)  ij  , 

vkykm;; - ^ - =  Km),  meft.  (2.16) 

ij 

Integrating  (2.15)  over  ft,  we  obtain 


Jv(y(m))  I  J(m)  I  dcf  =  Jl(m)d<7. 


Applying  to  the  Integral  on  the  left  the  known  formula  for  change  of 
variables,  we  obtain 

J  V(y)|  dt^l  =  |l(m)d<?  (2.17) 


This  formula  expresses  the  energy  conservation  property  of  the  reflector 
system.  It  is  known  in  the  literature  as  the  "balance"  equation;  see  [5], 
[11. 


3.  The  reciprocal  reflector 

3.1.  The  representation  of  F  as  a  graph  of  a  function  p  over  the 
"Input"  domain  ft  allows  the  construction  of  another  reflecting  surface 

F*  naturally  associated  with  F.  The  surface  F*  is  constructed  so  that  ft 
becomes  the  far-fleld  domain  while  w  becomes  the  "Input"  domain. 

Define  the  surface  F*  by  the  map 
-  r*  =  Vp  +  (p-'Ti)m,  m  £  ft, 


(3.1) 
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p  =  1/p,  -q  =  (p2  +  I  Vp|^)/(2p).  (3.2) 

It  Is  shown  in  [15]  that  r*  has  the  following  properties: 

(i)  q*  :=  5r*/9u’  =  q^jej’^Crrij^  -  (pj^/p)m),  i  =  l,2, 
where  q^j  =  V^jp  +  (p-iq)e|j; 

(ii)  the  vector  field 


^p  +  pm 


satisfies  the  relation  <rj  ,N>  =  0,  and  therefore,  if  r  is  an 

immersion,  then  N  is  the  unit  normal  vector  field  on  F*; 

(iii)  put  Q  ->  S^,  ^(m)  =  r*(m)/ri(m).  Then  obviously  N,  and 
y  are  coplanar  and  -<^,N>  =  <m,N>,  that  is,  the  law  of  reflection  is 

satisfied.  However,  this  time  the  reflector  F*  is  parametrized  by  points 
in  H. 

Explicit  computations  give: 


-r 


(3.3) 


(3.4) 


3.2.  Proposition.  Suppose  the  map  ^  defined  by  F*  is  a 
diffeomorphism  of  sT  onto  some  c  s.  Assume  that  the  map  7$:  Q  —> 
is  a  diffeomorphism.  Then,  t?(m)  =  ^(m)  for  all  m  e  Q,  and, 
consequently,  Q'  =  o>. 
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Proof.  Since  m  is  a  reflected  direction  relative  to  F  ,  we  have  by 
Snell's  law 


m  =  ^  -  2<$,N>N. 


(3.6) 


We  compute,  using  (3.3)  and  (3.5), 


iq<^,N>  =  <r*,N>  = 


I’pI 


2  .2  ,  ift.  .2  „.2  .2  .  i«.  |2 


p  n/P  »  |Vp 


2PWP^  *  |Vp| 


Also,  'r\<^,N>N  =  (Vp  -pm)(2p)"^  and 

Tim  =  Ti^  -  (Vp  -  pm)  p'2,  (3.7) 

Since  y  =  m  -  2  <m,n>n,  we  obtain  from  (1.1) 

Tig  =  ir\m  -  (pm  -  Vp)p"^,  (3.8) 

Then  from  (3.7),  (3.8)  we  get 

y(m)  =  ^(m)  for  all  m  €  (3.9) 

The  proposition  is  proved. 

3.3.  It  follows  from  (3.9)  that  =  id:  w  ->  w.  Therefore, 

=  J(^). 

This  is  also  confirmed  by  a  direct  computation.  Namely,  it  was  shown 


in  115]  that 


1? 


JCa  =  M(p)  = 


detTv.  D  +  (p-in)8.. 

ij  ij 

'q^detCe..) 


Then 


M(1/p) 


det 

1 

1 

2pV  p  -  ■ 
^  if 

(p2  -  |Vp|^)e..-  4p.pj  • 

(p^  *  Ivpl^) 

det(e  ) 
i] 

Comparing  with  (2.13),  we  see  that 

M(1/p)  =  G(p)  (3.10) 

This  formula  will  be  useful  in  several  instances. 

The  surface  F*  is  called  the  reciprocal  reflector.  With  the  use  of 
complex  structure  on  S,  it  has  been  described  in  [2]  and  {Wl,  chapter  2, 
section  6. 


4.  Differential  equations  of  the  problem 

4.1.  First  Order  System.  In  this  approach,  followed  by  Galindo- 
Israel  et  al.  [4-6]  the  analytic  formulation  of  the  problem  is  based  on 
the  following  considerations. 

The  equation  (1.8)  written  in  component  form  is  a  system  of  three 
partial  differential  equations  of  the  first  order.  However,  since  m  and 
y  are  unit  vectors,  we  may  reduce  (1.8)  to  a  system  of  two  equations  as 
follows. 

Observe  that  by  (1.8) 
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<y,m>  =1-2 


|Vp|2  +  p2  * 


(4.1) 


since  p  >  0,  we  may  put  v  :=  Inp.  Then  V|^  :=  dv/0u*^  =  p|^/p  and  we  get 


<y,m>  =  1  - 


2  I Vv|2 -  1 


+  Vv  2  1  +  yy  2 


(4.2) 


Again  from  (1.8)  we  get 


PPi 


2v. 


<y,m  >  =  2—5 — - 

^  k  l7p|2  +  p2 


1  +  |Vv|2 


,  k  =  1,2. 


(4.3) 


We  use  (4.2)  to  solve  for  |  Vv  1 2  +  i.  The  result  Is 


t^v|2  +  1  = 


1  -  <y,m> 


Then  we  obtain  from  (4.3) 


V 


k 


<y,m  > 

1 i  >< 

1  -  <y,m>’ 


1,2. 


(4.4) 


Thus,  If  both  domains  52  and  w  are  given  and  we  are  given  the  vector 
field  y(m),  then  by  solving  (4.4)  for  v  we  can  recover  p  and, 
therefore,  the  reflector  surface  F.  However,  for  the  system  (4.4)  to  be 
solvable,  an  Integrablllty  condition  must  be  satisfied  13],  Ch.  I.  In  this 
case,  the  condition  Is 

V 1 2  =  V2 1  (4.5) 

where  v^|^  =  92y/9ul9uk^  _  ^^2. 

Therefore,  for  a  reflector  surface  to  exist,  the  following  condition, 
derived  from  (4.4),  must  be  satisfied. 
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(<y2,m2>  -  <y2,mi>)(1  -  <y,m>)  = 

<y,m2><yi,nn>  -  <y,mi><y2,m>.  (4.6) 

In  a  different  form,  this  condition  Is  given  In  [7]  and  In  14-6]. 

Following  [5-6],  one  may  formulate  the  problem  of  synthesizing  a 
reflector  as  a  question  oj^ solvability  of  the  system  (4.4)  supplemented 
with  equation  (2.17)  with  prescribed  Q,  w,  and  functions  I:  Q  ->  (O,oo) 
and  V:  ->  (O,oo).  Since  In  terms  of  the  function  v  we  have  from  (1.8) 

y(m)  =  m  -  2  — (4.7) 

1  +  Vv  2 


the  system  (4.4)  and  equation  (2.16)  appear  as  a  first  order  system  of 
PDF's.  However,  as  the  formula  (2.13)  shows,  the  second  derivatives  are 
Involved  In  the  expression  for  the  Jacobian. 

In  a  series  of  papers,  Galindo-Israel  and  his  coauthors  show  that 
under  some  special  circumstances,  one  can  use  heuristic  arguments  for 
constructing  approximate  solutions  of  this  first  order  system  .  As 
mentioned  earlier,  their  results  rely  on  numerical  calculations. 

4.2.  In  some  special  circumstances,  the  Integrablllty  conditions 
(4.6)  can  be  easily  verified.  For  example,  consider  the  case  when  Q  Is  a 
circular  domain  with  center  at  the  North  pole  of  S  and  to  Is  also  a 
circular  domain  with  the  center  at  the  South  pole.  Let  the  vertical  axis 
Ox  pass  through  the  North  pole  and  (k,^)  be  spherical  coordinates  on  S 
such  that  0<c<<TT,  0<^<  2tt,  with  ^  =  0  corresponding  to  the 
positive  direction  of  the  z  axis. 

Suppose  that  the  map  7$  Is  such  that  y(m)  describes  w  as  a  surface 
of  revolution  about  the  z  axis  and  y(m)  Is  obtained  from  m  by  rotating 
m  In  the  plane  passing  through  Oz  and  m;  scs  the  fig.  below. 
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A  Z 


(x,y) 


We  want  to  check  that  (4.6)  is  satisfied.  Since  m^:=  0nn/9^  is 

perpendicular  to  the  plane  spanned  by  m  and  the  Oz  axis,  we  have 
<y(m),m^>  =  0.  Because  Q  and  (o  are  rotationally  symmetric  about  the 

Oz  axis,  <y(m),m>  =  f(^)  for  some  function  f.  Then  <y^(m),m>  =  - 

<y(m),m^>  =  0  and  substituting  into  the  right  hand  side  of  (4.6)  we 

obtain  (i  =i<-->c<,  k=2<-->^) 

<y,m^><y^,m>  -  <y,m^><y^,m>  =  0 
Differentiating  <y,m^>  in  a,  we  obtain 

But  it  is  easy  to  check  that  is  perpendicular  to  m  and  to  the  Oz 

axis.  Hence,  it  is  perpendicular  to  y(m).  Then  <yj^,m^>  =  0. 

Similarly,  one  checks  that  <y^,mp^>  =  0.  Thus  (4.6)  is  satisfied. 

In  general,  one  does  not  know  3  priori  the  vector  function  y(m),  but 
if  such  a  function  exists,  it  must  satisfy  (4.6). 

4.3.  The  Monge-Ampere  Equation.  In  this  approach,  one  attempts 
to  solve  the  equation  (2.16).  let  us  formulate  the  problem  precisely. 
Suppose  we  are  given  two  domains  Q  and  o;  on  S  and  two  positive 
functions  I:  Q  — >  (0,oo)  and  V:  w  ->  (0,oo).  The  problem  consists  In 
finding  a  solution  p  >  0  of  the  equation 
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V(y(m))(G(p))(m)  =  I(m),  m  e  Q,  (4.8) 

or 

V(y(m))(G(p))(m)  =  -I(m),  m  c  Q,  (4.8)' 

subject  to  the  boundary  condition 

dQ  ->  (4.9) 

In  (4.8)  y(m)  is  expressed  in  terms  of  p  as  in  (1.8).  The  boundary 
condition  (4.9)  is  understood  as  the  requirement  that  the  boundary  is 
mapped  homeomorphically  onto  do;,  but  the  map  is  not  specified  point- 
wise.  One  can  show  that  if  7$  is  specified  point-wise  on  dQ,  then  the 
problem  is  overdetermined. 

The  data  I,  V,  and  o;  cannot  be  arbitrary  since  the  energy 
conservation  requirement  (2.18)  must  be  fulfilled.  The  latter  can  be 
rewritten  as 

Jv(y)da(y)  =  Jl(m)da(m)  (4.10) 

«  n 

The  operator  G(p)  in  (4.8)  and  (4.8)'  is  of,  so  called,  Monge-Ampere 
type  and  whether  it  is  elliptic  or  hyperbolic  depends  on  the  class  of 
functions  on  which  it  is  considered.  More  precisely,  G  will  be  positively 
(negatively)  elliptic  on  any  p  for  which  the  matrix 

(2p^,jP  -  (p2  -|Vpr)e,j  -  4piPjl  Maij(p)l 

is  positive  (negative)  definite.  G  will  be  hyperbolic  on  such  p  for  which 
lajj(p)l  is  indefinite  but  nondegenerate.  Since  la^j(p))  is  a  2  x  2  matrix, 

G  is  elliptic  if  and  only  if  detla^jfp)]  >  0  in  and  hyperbolic  if 

det[ajj(p)I  <  0  in  Q.  Respectively,  equation  (4.8)  corresponds  to  elliptic 

solutions,  both  positive  elliptic  and  negative  elliptic,  and  additional 
restrictions  on  p  needs  to  be  imposed  in  order  to  specify  one  of  them. 
The  following  simple  example  illustrates  this  situation.  Take  a 
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surface  F  given  by  the  function  p  -  c  =const.  >  0.  Then  F  is  a  sphere  of 
radius  c,  a^j(c)  =  G(p)  =  1  (see  (2.13))  and  the  matrix  [a^j(c)]  is 

negative  definite. 

Consider  now  a  surface  F  which  is  plane  tangent  to  S  at  the  North 
Pole.  One  can  calculate  p  and  similarly  compute  the  corresponding 
[a|j(p)].  However,  a  simpler  way  is  to  observe  from  (2.12)  and  (2.13) 

that  la^j(p)](p^  +  |vp|^)"^=  [2b^j<m,n>  +  pe^jlp”^  and  then  recall  that 

for  a  plane  b|j  =  0  [9].  Therefore,  in  this  case  la^j(p)]  is  positive 

definite,  since  [e^j]  is  positive  definite. 

Finally,  note  that  the  equation  (4.8)'  corresponds  to  those  p  on 
which  G  is  hyperbolic. 

Thus,  prion,  without  considering  G  on  a  particular  p,  one  cannot 
describe  the  type  of  G  even  though  I/V  is  positive.  We  emphasize  this 
standard  point  only  because  in  [16]  it  was  claimed  that  the 
corresponding  equation  is  elliptic. 

Let  us  now  show  that  on  elliptic  ‘solutions”  the  Jacobian  J(m)  is 
always  positive.  Indeed,  if  the  map  ->  o)  preserves  the  orientation, 
the  products  <m,mjX  m2>  and  •<y,U|Xy2>  have  the  same  sign.  Since  by 

our  assumption  <m,m|Xm2>  >  0  in  the  product  <y,yixy2^  >  0  in 

and  consequently,  J(m)  >  0  in  In  view  of  (1.5)  and  (2.13),  J(m)  and 
(G(p))(m)  must  have  the  same  sign.  Therefore,  (G(p))(m)  >  0,  which 
means  that  J(m)  is  positive  only  on  such  p  for  v/hich  G  is  elliptic. 

Similarly,  one  shows  that  J(m)  <  0  on  those  p  on  which  G  is 
hyperbolic. 


5.  Radially  symmetric  case 

5.1.  Now  we  consider  the  special  case  in  which  the  reflector  has 
axial  symmetry.  In  this  case,  the  PDE  (4.8)  reduces  to  an  ordinary 
differential  equation  for  which  we  can  find  explicit  solvability 
conditions.  It  will  be  shown  below  that  the  situation  here  is  very 
similar  to  the  case  considered  in  [15],  and  we  follow  this  work  closely. 

It  will  be  convenient  in  this  case  to  use  spherical  coordinates  a, 
where  -tt/2  <  k  <  tt/2,  0  <  ^  <  2tt.  Assume  that  both  domains  Q  and  to 
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are  circular  with  centers  being  respectively  the  North  and  South  Poles. 

I  oc  <  K  <  TT/2},  cx  €  (0,tT/2)  , 

=  {  (cx,^)  I  -tT/2  <  o<  <  o<},  (X  €  (-Tt/2,0). 

In  these  coordinates, 


0 

cos^cx 


We  are  looking  for  elliptic  solutions  of  (4.8)  in  the  form  p  =  p(cx). 

We  introduce  a  new  unknown  function  p  =  1/p,  assuming,  as  always, 
that  p  >  0.  As  it  was  shown  in  section  3,  the  operator  G(p)  =G(1/p)  = 
M(p).  Using  the  expression  for  M(p)  in  the  axially  symmetric  case  as  in 
[151,  section  3,  we  obtain  from  (4.8) 

11/  /  \nm/  n  li/  /  ^2pp  +  p^  -  p^)(-2pptan  c<  +  p^  -  p) 

V(y(m))M(p)  s  V(y(m))  — — — - - - !— =  I(m),  (5.1) 

(p2  +  p2)2 

oT  <  cx  <  tT/2, 


where  p  =  dp/dcx  p  =  d^p/dcx^. 

In  order  for  p(cx)  to  be  smooth  in  Q,  we  need  to  impose  the  condition 
p(Tt/2)  =  0  (5.2) 

which  is  equivalent  to 

P(tt/2)  =  0.  (5.3) 

We  now  set  up  the  boundary  conditions.  This  amounts  to  describing 
(4.9)  in  analytic  form.  If  one  prescribes  1$:  dO.  ->0a>  point-wise,  then 
the  problem  is  overdetermined  (  cf.  [5],  [6],  [7],  [8]).  A  way  to  relax  this 
restriction  is  to  require  that  maps  9^2  ohto  0w  homeomorphically,  but 
without  specifying  the  map  point-wise. 
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Thus,  if  y  =  m  €  then 

<m,(0,0,1)>  =  cos(Tr/2  -  a)  =  sin^  on  dQ, 

<y,(00,-l)>  =  cos(tt/2  -  a)  =  -sinoc  on  du). 

On  the  other  hand,  by  (1.8)  we  have,  taking  into  account  that|Vp|2  =  p2, 

2p(pm  +  Vp) 
y  =  m - rr - 5 - . 

p2  +  p2 

Then 

-  sinS  =  -  sin5-  +  ..  . ..  sini^  -  ,  <^p,(0,0,-  1)>  . 

p^  +  p'^  +  p^ 


Since  Vp  =  p  9m/9o<,  we  have  <Vp,(0,0,-1)>  =  -pcosc<  and 


.  „  _  2pp 

-  sincx  =  — - - —  sinc<  +  - - -cosa 

p2  .  p2  p2  .  p2 

or. 


p2(sino<  -  sinoc)  =  p2(sino<  +  sinoc)  +  2pp  cosoc  =  0  (5.4) 

where  p  and  p  are  evaluated  at  cc  =  a.  This  boundary  condition  is 
identical  with  (2.3)  in  [15].  In  the  special  case  when  the  illumination 
pattern  is  required  to  be  uniform,  that  is,  V(y(m))  =  const  =  Vq,  then  the 
problem  (5.1),  (5.3),  and  (5.4)  is  analytically  the  same  as  the  one  in  [15] 
and  we  can  formulate  the  following  result  (cf.  Theorem  3.5  in  [15]). 

Theorem  5.2.  Let  I(c<)  be  positive  and  continuous  on  [O.tt/2].  Let  oc 
be  anu  number  in  the  interval  (-tt/2,0)  and  k  the  solution  of  the  equation 


2 

1  +  sinoc  =  J  I(r)costrdr,  cr  e  (0,Tt/2) . 
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Then  for  each  choice  of  p(o(),  p(c<)  such  that 
p(cr)  ■  2  ■  2 


or 


tT  SL  +  ^ 
p(or)  “  2  ■  2  2 

there  exists  a  unique  solution  p  >  0  of  (5.1),  (5.3),  (5.4)  of  class 
CUc<,TT/2]nC^  (c<,TT/2]. 

5.3.  When  the  function  p  is  constructed,  then  we  return  to  the 
function  p  =  1/p  and  consider  r(m)  =  p(o«:,^)m(c<,^)  where  we  set 
p(o<,^)  s  pioc).  The  vector  function  r  defines  the  reflector  surface. 
Clearly,  the  map  ^  is  defined  and  J(m)  >  0.  Thus  is  a  local 
diffeomorphism. 

6.  An  algorithm  for  computing  the  output  power  density. 

It  follows  from  formulas  (2.16)  and  (2.12)  that  the  computation  of 
the  output  power  intensity  reduces  to  computation  of  the  Jacobian  of 
the  map  '6.  In  this  section  we  provide  a  discretization  of  (2.12)  and 
several  computational  examples.  We  preserve  here  the  notation  from 
previous  sections. 


ALGORITHM  FOR  THE  COMPUTATION  OF  G(p) 
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Ti  (a+sh  Ik)  =  -  r(a+(s+l)h./k) 

2h 

T,(S«h,nc)  = 

c.  if  s  =  0  (i.e.  in(a+sh,/k)  lies  on  the  boundary  of  Q)  we  make  a  one-sided  approximation 
of  Si(a-Hsh,flc)  and  Tj  (a-t-sh,/k): 

S|(atsli.lk) .  -  ni(atsh.ft)  ^ 


s,(5+sh,/k)=. 
T,(5+sh,/k)  = 


T2(aH-sh,/k)  = 


r(a-Hsh./-Hl)  -  r(a-Hsh,/-l) 
2k 


h 


Remark:  T,  (a,/k)  best  approximates  a  tangent  vector  at  the  r(a-(-— ,/k). 


2.  For  0  Ss  <  s-1  and  1  <  /  <7-1  estimate  n(a-t-sh,/k)  as 


Tj  (a+sh,/k)xT2(a-t-sh,/k) 


ITi(a+sh,/k)xT2(a+sh,/k)| 

3.  We  approximate  the  directional  derivatives  of  n(a+sh,flc)  with  respect  to  Ti(a-HSh,/k)  and 
T2(a+sh,/k)  once  again  in  3  ways: 

a.  if  s  =  s-1 


_  n(a+(s-l)h,(Y+/  mod(O)lc)  -  n(a-i-(s-l)h,/k) 

n,  (a+sh,/k)  - - — - , 

n(a+(s-l)h,(  j-t-/  mod(0)k)  -  n(a-i-(s-l)h,(-j-+/  mod(0)k) 
n2(a+sh,/k)  = - — - . 

b.  if  1  <  s  <  s-1 
n,(5+sh.nc)  = 

2h 

c.  if  s=  1,  we  use  the  remark  after  step  Ic  to  conclude  that  since  Ti(a,/k)  best  approxi¬ 
mates  the  tangent  vector  at  rCa+y./k),  n(a,/k)  best  approximates  the  normal  at 
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r(a+|-,/k): 

n  ^  n(5+(s-l)h,/k)-  n(a-f(sH-l)h.yk)^ 

l.jh 

.,(5+sh.ik)  = 

4.  Approximate  the  metric  at  m(a+sh,/k)  with  respect  to  Si(a+sh,/k)  and  S2(a+sh,/k)  as 
eij(a+sh,/k)  =  <Si(a+sh,/k),Sj(a+sh,/k)>  for  i  =  1,2  and  j  =  1,2. 

5.  Approximate  the  second  fundamental  form  at  r(a+sh,flc)  as 


bij(a+sh,/k)  =  -<ni(a+sh,/k),Tj(a+sh,/k)>fori  =  l,2andj=  1,2. 


Evaluate  G(p(a+sh,/k))  for  0<s<s-2  and  0  </</-! 


det(eij  (a+sh,  /k))  ‘  det  2 


,  bii  (a+sh,  /k)<m(a+sh,/k),n(a+sh,flc)> 

O _ i _ 


p(a+sh,flc) 


+eij(a+sh,/k)^ . 


We  used  this  algorithm  for  examples  where  G(p)  is  known.  These  were  for  p(a,P)  =  1  (the  unit 

2 

sphere  with  the  center  at  the  origin),  for  p(a,P)  =  csc(a)  (the  plane  z  =  1)  and  for  p(a,B)  =  - — :  .  . 

l+sm(a) 

(the  paraboloid  with  equation  z  =  1  -  -^(x^  +  y^)  with  its  focus  at  the  origin). 

For  all  three  surfaces  we  took  a  =  -^ ,  s  =  25  andT  =  24. 

4 

Example  1:  For  the  unit  sphere  G(p)  s  1.  Our  algorithm  had  a  maximum  error  of  8.9x10“^, 
which  occured  for  s  =  1. 

Example  2:  For  the  plane  G(p)  s  1.  Our  algorithm  computed  this  with  no  error  since  it  com¬ 
puted  the  second  fundamental  form  as  precisely  0. 

Example  3:  For  the  paraboloid  G(p)  s  0.  Our  algorithm  had  a  maximum  error  of  7.6x10"’, 
which  occured  for  s  =  1. 

We  have  also  used  this  algorithm  on  offset  examples. 
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On  one  direct  problem  in  the  reflector  antenna  theory 


by 

Vladimir  Oliker’ 
and 
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Abstract 

This  paper  is  the  second  in  a  series  of  papers  in  which 
we  use  differential  geometric  methods  to  investigate 
systematically  the  problem  of  synthesis  of  single  and  dual 
reflector  antennas.  In  our  first  paper  we  considered  the 
direct  problem  for  a  single  reflector  and,  in  particular, 
established  rigorously  existence  of  radially  symmetric 
solutions  in  the  case  when  the  data  is  radially  symmetric.  In 
present  article  we  prove  existence  of  reflectors  solving  the 
direct  problem  in  the  case  when  the  data  is  not  radially 
symmetric  but  close  to  such  in  some  Holder  norm. 


1.  Introduction 


This  paper  is  the  second  in  a  series  of  papers  in  which  we  study  the 
theory  and  numerical  methods  in  synthesis  of  reflector  antennas.  We  use 
the  geometric  optics  (GO)  approximation  to  describe  and  study  the 
problem.  A  brief  history  of  the  work  of  other  researchers'  in 
constructing  reflector  antennas  using  GO  can  be  found  in  our  paper 
(ONPl. 

'  Research  sponsored  by  the  AFOSR  (AFSC)  under  Contract  F49620-91-C-0001.  The  United  States 
Government  is  authorized  to  reproduce  and  distribute  reprints  for  government  purposes  notwithstanding  any 
copyright  notation  hereon. 

’  The  author  is  a  graduate  student. 
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Building  upon  our  previous  work  [ONP]  and  10],  we  show  now 
existence  of  nonradially  symmetric  solutions  do  direct  problem  in  the 
case  where  the  density  of  distribution  of  reflected  rays  is  uniform  and 
the  density  of  the  distribution  of  the  incidence  rays  is  not  radially 
symmetric  but  close  in  a  certain  norm  to  a  radially  symmetric  function. 


2.  Preliminaries  and  the  main  result. 

2.1  We  begin  by  recalling  the  notation  used  in  [ONPj.  In  three 

dimensional  space  let  S  be  a  unit  sphere  centered  at  the  origin  0  of  a 
Cartesian  coordinate  system.  Let  be  a  domain  on  S,  Q  its  closure,  and 
m  a  unit  vector  with  endpoint  in  Q.  Let  p  be  a  positive  function  Of  class 

C^(Q)ncHQ)  and  set 


r(m)  =  p(m)m.  (1) 

The  map  r  defines  a  surface  F  projecting  radially  from  0  univalently 
onto  Q.  Denote  by  n  the  unit  normal  vector  field  on  F  and  assume  that  F 
is  oriented  so  that  <m,n>  >  0  everywhere  on  F.  Here,  <  ,  )  denotes  the 

usual  scalar  product  in  R^. 

If  a  light  ray  originates  at  0  in  the  direction  m,  reaches  F,  (the 
reflector)  and  is  reflected,  then  we  can  put  in  correspondence  with  m  a 
unit  vector  y  e  s  parallel  to  the  reflected  ray.  Thus,  we  have  the 
"reflector"  map  ?2">S,  y  =  l^(m).  The  laws  of  geometric  optics  tell  us 
that  the  vectors  m,  n(m),  and  tf(m)  lie  in  one  plane  and  ?f(m)  =  m  - 
2<m,n)n.  The  image  Z>  =  is  called  the  "far  field". 

Let  u  =  (u^u^)  be  some  smooth  local  coordinates  on  S  such  that  Q 
lies  inside  one  coordinate  patch  (it  is  always  assumed  that  Q  ^  S).  Then 

m(u)  s  m(u^u^)  is  a  smooth  vector  valued  function  giving  the  position 
vector  of  any  point  in  Q.  For  that  reason  m(u)  is  viewed  as  a  unit  vector 
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in  and  also  as  a  point  in  H.  As  usual,  we  put  f(m(u))  s  f(u)  for  any 
function  f:  S— >R.  Everywhere  in  the  paper  we  use  the  range  of  indices 
1<i,j,k,..<2. 


The  first  fundamental  form  e  =  e^jdu^duJ  of  S  has  coefficients 

e^j  =  where  m|  =  9m/du\  Here  and  everywhere  below  the 

convention  about  summation  over  repeated  lower  and  upper  indices  is  in 
effect.  The  matrix  [e|j]  is  symmetric  and  invertible;  its  inverse  is 

denoted  by  [e’j]. 

Set 

Vp  =  p^e^Jmj. 

Then 

|Vp|2  =  <9p,Vp>  =  p,pje*j. 


It  is  shown  In  lONPl  that  for  any  p  «  p  >  0  In  the  vector  field 


pm  -  Vp 


(2) 


is  the  unit  normal  vector  field  on  the  surface  F  defined  by  (1).  In 
addition,  <r,n)  >  0  on  F.  Further,  if  y  =  ?f(m)  is  the  unit  vector  parallel  to 
the  reflected  ray  corresponding  to  m  then 


y  =  m  - 


p(pm  -  Vp) 

I  2 

p  *  I’^pI 


(3) 


2.2  In  the  difiect  r^ffector  p/vd/e/n  the  following  data  is  given:  the 
domains  Q,  w  on  S,  a  positive  function  I:  — >(0,oo),  and  a  positive 
constant  Vq.  We  have  to  find  a  reflector  F  subject  to  the  requirements: 
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(1)  the  rays  originating  at  0  and  going  through  points  of  Q  project 
F  univalently  onto  Q; 

(ii)  <i>  is  the  far  field; 

(iii)  for  the  given  density  Km)  of  the  incidence  rays  the  density  of 
the  distribution  of  rays  reflected  of  F  in  the  direction  ?^(m)  is  uniform 
and  equal  to  Vq; 

(iv)  finally,  it  is  also  natural  to  require,  in  this  setting,  that  F  be 
such  that  the  reflector  map  is  a  diffeomorphism  of  Q  onto  (o. 

It  is  shown  in  [ONP]  that  in  terms  of  function  p  the  Jacobian  of 
the  map  7$  is  given  by 


G(p)  := 


Vp|  ) 


det  [2pV  p  -  (p^  -  I  Vp|  )e..  -  4p.p.] 

r  T  I  r  I 

det  [e.  .1 


(4) 


Thus,  in  order  to  find  a  reflector  F  satisfying  (iii)  we  need  to  solve  the 
PDE  (cf.  (ONPl,  subsection  2.5) 

G(p(m))  =  I(m)/VQ,  in  (5) 

If  the  reflector  map  7$,  determined  by  F,  saUsfies  (i)  and  (iv)  then  it  is 
necessary  that 


7i:  dQ—>d<i>.  (6) 

By  (2),  7i  is  expressed  in  terms  of  p,  and,  therefore,  (6)  is  a  condition  on 
p  on  Hence,  we  may  view  it  as  a  boundary  condition  to  be  satisfied 
by  solutions  of  (5).  However,  if  (6)  is  understood  as  a  pointwise 
condition,  then  the  problem  (5),  (6)  is  overdetermined  and,  in  general, 
will  not  have  a  solution.  For  this  reason  (6)  is  treated  as  a  requirement 
that  the  boundary  of  be  mapped  onto  the  boundary  of  <s>. 
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2.3  In  order  to  formulate  our  main  result  we  need  some  more  notation. 


Let  ((x,^),  -Tt/2  <  (X  <  tt/2,  0  <  ^  <  2Tt,  be  the  spherical  coordinates 
on  S.  Let 

=  {((x,^)  I  (X  <  (X  <'rT/2},  K  G  (0,tT/2), 
w  =  {(cx,^)  I  -tT/2  <  oc  <  «},  «  G  (-Tr/2,  0). 

Let  r  be  a  positive  function  of  class  cMo,tt/2].  Put 

E(c<  ,(x  )  =  f  r(r)cos  X  dr. 

12  ^ 

I  ^ 

For  5  G  (0,1)  put 

O.S  _  .  , 

H  =  {  I  G  C  (^^)  I  |l  dcj  =  2TtV  (1  +  sin  «)},  (7) 

where  da  is  the  area  element  of  S. 

Theorem  A.  Suppose  that  Q,  co  and  I  are  such  that 

E(k,  tt/2)  =  Vq(1  +  sin  a).  (8) 

Then  there  exists  an  €  >  0  such  that  for  any  I  g  H  satisfying  ill  -  I  ii  <  €, 
where  li  ii  denotes  the  Holder  norm  in  C^»^($2),  the  equation 

G(p)  =  I/Vq  in  (9) 

with  the  boundary  condition 

(?J(m),(O,O,-1)>|0j2  =  geodesic  radius  of  To 
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(10) 


admits  two  classes  of  solutions  p  e  Within  each  class  the 

solution  is  determined  uniquely  up  to  a  positive  multiplicative  constant. 
Furthermore,  the  corresponding  reflector  surfaces  satisfy  the  conditions 
(i)-(iv)  as  described  previously. 

Remarks. 

(a)  The  condition  I  e  H  together  with  (7)  expresses  the  requirement 
that  the  data  must  satisfy  the  energy  conservation  law  (cf.  [ONP], 
subsection  2.3). 

(b)  If  for  (x  <  cx  <  tt/2,  0  <  ^  <  2tt,  we  set  Kk,^)  =  Koc)  then  (8) 
implies  that  Fe  H.  It  is  shown  in  [ONP],  Theorem  5.1,  that  for  Q,  (x>,  F, 
and  Vq  satisfying  (8)  one  can  construct  two  classes  of  radially 

symmetric  solutions  to  (9),  (10)  (with  F  on  the  right  hand  side  of  (9)) 
and  in  each  class  two  solutions  differ  at  most  by  a  positive 
multiplicative  constant.  The  role  of  the  function  F  and  any  of  the 
corresponding  radially  symmetric  solutions  is  that  these  special 
solutions  serve  as  approximations  to  nonradially  symmetric  solutions  of 
(9),  (10).  In  this  article  this  connection  remains  in  the  background  and  is 
not  explicitly  shown,  since  the  proof  of  Theorem  A  is  based  on  a 
reduction  of  the  direct  problem  to  the  inverse  problem  investigated  in 
[Oj.  For  the  inverse  problem  this  connection  is  explained  in  detail  in  [0]. 

(c)  Let  I  6  HnCkfi).  Put 

Im(cx)  =  ^  I  I(oc,^)d^ 

0 

It  is  easy  to  see  that  the  function  1^^  satisfies  all  the  conditions 

imposed  on  I.  Thus,  the  hypothesis  regarding  existence  of  the  function  F 
can  be  replaced  by  the  requirement  that  ill  -  1^11  <  e. 

3.  Proof  of  Theorem  A 

Step  1.  First  we  show  that  there  exist  two  distinct  classes  of 
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functions  satisfying  (9),  (10).  Following  lONP],  section  3,  we  introduce  a 
new  unknown  function  p  =  1/p.  Then  (9)  becomes 


M(p) 


det  IV. .p  +  (p-in)e. .1 
U  ij 

det  [e  ] 
ij 


in  Q, 


(11) 


where  =  (p^  +  |Vp|^)/2p.  Since  (10)  is  not  written  explicitly  in 
terms  of  p,  it  will  remain  the  same,  but  in  order  to  emphasize  that  we 
are  using  the  function  p,  we  rewrite  it  as 


<?fp(m),(0,0,- 0)  =  geodesic  radius  of  <x>.  (12) 

With  any  p  >  0,  which  is  a  solution  of  (11),  (12),  we  associate  the 
reciprocal  reflector  F*,  that  is,  the  surface 

^  ^  _ 

-r  =  Vp  +  (p-'T\)m,  m  e 


We  check  by  a  direct  computation  that  |r*|  =  T[  and  ^  ;=  r*/'T\  =  ^6^  in 

a.  Thus,  the  problem  (11),  (12)  is  identical  to  the  problem  (3.1),  (3.2)  in 
[0]  (with  the  roles  of  Q  and  <jl)  interchanged).  On  the  other  hand,  if  in 
Theorems  2.1  and  3.1  in  [01  we  set  T=  !?Vq  and  interchange  the  roles  of 

a  and  oc  then  the  hypotheses  of  our  Theorem  A  imply  the  hypotheses  of 
Theorems  2.1  and  3.1  in  [0].  Consequently,  we  conclude  that  the  problem 

(11),  (12)  admits  two  classes  of  positive  solutions  in  C^(?2)  and  within 
each  class  the  solution  is  unique  up  to  a  positive  multiplicative 
constant.  Furthermore,  in  case  I  =  F  the  solutions  of  (11),  (12)  are 
radially  symmetric. 

We  return  now  to  the  original  unknown  function  p  =  1/p  and  that 
completes  the  step  1  of  the  proof. 

Step  2.  Let  p  be  a  solution  of  (9),  (10)  and  F  the  surface  defined 
by  (1).  On  this  step  we  check  the  conditions  (i)-(iv)  in  subsection  2.2. 


7 


It  follows  from  the  discussion  in  sections  2.1  and  2.2  that  F  is  a 
reflector  surface  with  the  normal  field  given  by  (2)  and  the  far  field 
described  by  vectors  y  =  ?^(m)  given  by  (3).  The  condition  (i)  is  obviously 
satisfied.  Condition  (iii)  is  satisfied  because  (9)  is  satisfied.  It  remains 
to  check  (ii)  and  (iv). 

Since  the  Jacobian  J(2^)  =  G  =  I/Vq  >  0  in  Q,  the  map  is  a  local 

diffeomorphism.  Because  of  (10)  ^  maps  Q  into  a>.  By  Corrolary  4.7  in 
[KN],  since  H  is  compact,  the  map  is  a  covering  projection.  On  the 
other  hand,  if  is  not  a  global  diffeomorphism  of  Q  onto  a>  then,  since 
J(?J)  >  0, 

area  of  <jl)  =  2tt(1  +  sin  &)  <  |j(tJ)d<?  =  jG(p)dcf  =  Jl  dcJ. 

The  latter  contradicts  the  hypothesis  that  I  e  H.  Now  the  conditions  (ii) 
and  (iv)  are  also  verified  and  the  theorem  is  proved. 
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